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THE  VISIBILITY  OF  STATIONARY  AND  MOVING  TARGETS  IN  THE  PLANE 
SUBJECT  TO  A  POISSON  FIELD  OF  SHADOWING  ELEMENTS* 


by 

M.  Yadin  and  S.  Zacks 

Technicn,  Israel  Institute  o£  Technology  and 
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ABSTRACT 


A  methodology  for  an  analytical  derivation  of  visibility 
probabilities  of  n  stationary  target  points  in  the 
plane  is  developed  for  the  case  when  shadows  are  cast 
by  a  Poisson  randan  field  of  obscuring  elements.  In 
addition  formulae  for  the  moments  of  a  measure  of  the 
total  proportional  visibility  along  a  star-shaped  curve 
are  given.  The  general  methodology  is  illustrated  by 
some  examples  of  visibility  along  the  circunference  of 
half  a  circle,  when  the  obscuring  elements  are 
centered  within  a  concentric  annular  region.  The 
distribution  of  the  measure  of  proportional  visibility 
is  approximated  by  a  mixed-beta  distribution.  _ 


Key  Words:  V-Ui. bilAXy  probability,  stationary  and  moving  targets; 

measure.  o\  proportionaZ  oisibiJUty;  Poisson  Aandom  shadouiing 
process . 
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Research  Office. 


0.  Introduction 


U5e  present  paper  is  the  first  one  in  a  series  of  articles 
dealing  with  problems  of  visibility  of  targets  through  random 
fields  of  obscuring  elements  (trees,  bushes,  clouds,  etc.)  and 
similar  problems  of  detection  and  hitting  of  targets.  The 
problems  dealt  with  are  of  stochastic  nature.  The  exact  number, 
location  and  dimensions  of  the  obscuring  elements  are  unknown.  A 
probability  model  is  formulated  concerning  these  variables.  Given 
such  a  probability  model  it  is  required  to  determine 
probabilities  of  certain  events  and  distributions  of  certain 
random  variables.  To  illustrate  sene  of  the  problems  that  can  be 
solved  by  the  methodology  developed,  consider  the  following 
examples: 


Example  1:  Visibility  of  Stationary  Targets 

to  observer  is  placed  at  a  given  location  in  a  forest,  in 
order  to  detect  specified  targets  (vehicles,  animals,  etc.) .  Due 
to  the  random  location  of  the  trees  it  is  important  to  determine 
the  probabilities  that  individual  targets  are  observed  and  the 
distribution  of  the  number  of  targets  observed.  For  this  purpose 
one  has  to  determine  the  probabilities  that  any  specified  r 
points  out  of  n,  1  <  r  <  n,  are  simultaneously  visible.  The 
locations  of  the  targets  may  be  specified  car  random. 

Example  2:  Visibility  of  Moving  Targets 
A  target  is  moving  along  a  specified  path  C.  to  observer  is 
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located  at  a  point  0.  Trees  or  other  obscuring  elements  are 
distributed  between  0  and  C.  It  is  often  required  to  determine 
the  distribution  function  of  the  total  visible  portion  of  c. 


;  Detecting  and  Hitting  Moving  Targets 


A  hunter  is  located  at  a  point  0  in  the  forest  and  a  target 
is  moving  along  a  path  C ,  which  is  partially  obscured  by  the 
random  elements.  In  order  to  detect  a  target  the  hunter  should 
observe  it  uninterruptedly  for  tQ  units  of  time.  It  is 
interesting  to  determine  the  probability  of  detecting  a  target  at 
a  specified  location  or  anywhere  on  C.  In  addition,  after 
detecting  a  target  the  hunter  could  attempt  to  hit  it.  Each  such 
attempt  requires  t  units  of  time  in  which  the  target  is  visible. 
A  hitting  attempt  may  be  successful  with  probability  p.  If  an 
attempt  fails  the  hunter  can  attempt  again,  as  long  as  the  target 
is  visivle.  Once  the  target  disappears,  it  has  to  be  detected 
again.  The  total  number  of  possible  trials,  N,  depends  on  the 
length  of  the  visible  portions  of  C .  N  is  a  randan  variable,  and 
it  is  interesting  to  determine  the  distribution  of  N,  and  the 
probability  of  hitting  a  target. 


Example  4;  Biergy  Penetration  Through  Random  Fields 
A  light  beam  from  the  origin  0  is  focused  in  a  specified 
direction.  If  there  are  no  obscuring  elements,  a  given  sector  of 
the  plane  is  in  the  light.  On  the  other  hand,  a  portion  of  the 
sector  may  be  in  the  shadow  of  random  elements.  In  such  a  case 


part  of  the  energy  in  a  beam  does  not  penetrate  the  field.  The 
objective  is  to  determine  the  distribution  function  of  the  amount 
of  penetrating  energy. 

In  the  above  examples  we  restricted  attention  to  visibility 
problems  in  the  plane.  These  problems  can  be  readily  generalized 
to  visibility  in  three  dimensional  spaces,  where  the  obscuring 
elements  are  crowns  of  trees,  clouds,  sattelites,  stars,  etc. 
Problems  of  similar  nature  can  be  found  in  various  areas  of 
applied  science  and  technology  and  in  military  applications. 

The  literature  dealing  with  these  kinds  of  shadowing  prob¬ 
lems  contains  a  small  number  of  papers.  Chemoff  and  Daly  [2] 
discussed  the  problem  of  determining  the  distribution  of  the 
length  of  shadowed  and  of  visible  segments  an  a  straight  line  in 
the  plane.  Likhterov  and  Gurin  (6]  studied  the  probability  of 

detecting  a  moving  target  on  a  straight  line.  The  shadowing 
problem  can  be  reduced  to  a  coverage  problem.  This  reduction  is 
achieved  by  specifying  the  shadows  cast  by  random  elements  on  a 
target  curve.  In  a  previous  paper,  Yadin  and  Zacks  [12]  applied 
methods  of  coverage  probabilities  to  study  visibility  problems  on 
a  circle.  The  literature  on  coverage  problems  is  very  rich  (see 
Bckler  [3]).  The  papers  of  Ailam  [1],  Greenberg  [5],  Robbins 
[8,9]  and  Siegel  [10,11]  provide  methods  which  are  relevant  to 
the  study  of  shadewing  problems.  In  particular  the  methods 
developed  in  the  present  paper  and  those  of  [12]  are  based  on  the 
theory  given  by  Robbins  [9,  10] .  The  geometry  of  the  problem  may 
be  coopl ica ted,  and  a  reduction  to  a  coverage  problem  may  be 
inconvenient.  In  such  cases  it  is  often  easier  to  study  the 


geometry  of  the  shadowing  elements  rather  than  that  of  the 
shadows.  For  this  reason,  the  present  paper  deals  with  the 
shadowing  problems  directly. 


Die  fact  that  the  number  of  papers  on  shadowing  problems 
published  is  small  does  not  reflect  a  lack  of  interest  in  these 
problems  in  various  areas  of  application.  One  can  find  a  large 
number  of  technical  reports  which  treat  shadowing  problems  by 
simulations  or  by  fitting  simple  statistical  models  to  empirical 
data.  Such  "solutions”  do  not  died  light  on  the  intrinsic 
structure  of  the  stochastic  phenomenon.  Hie  methodology 
developed  in  the  present  paper  analyzes  the  stochastic  structure 
and  provides  methods  for  numerical  approximations  which  are  more 
accurate  and  require  less  computing  time  than  simulation 
procedures.  For  developing  the  general  methodology  we  construct 
in  Section  1  a  probability  model  of  random  shadowing  elements. 
Sections  2  and  3  deal  with  visibility  probabilities  of  single  and 
several  stationary  targets  (Example  1) .  Section  4  deals  with  a 
visibility  measure  for  moving  targets  (Example  2) .  The 

methodology  for  solving  problems  discussed  in  Example  3  and  4 
will  be  presented  in  forthcoming  papers.  The  methods  of  the 
present  paper  are  extended  in  [13]  to  solve  three  dimensional 
visibility  problems. 


1.  Random  Fields  of  Shadowing  Objects 


In  general,  a  shadowing  object  is  characterized  by  parameters 
of  location,  shape  and  size.  In  the  present  paper  we  restrict 
attention  to  disks  in  the  plane.  An  extension  to  spheres  in 
three  dimensional  spaces  is  given  in  [13] .  A  disk  in  the  plane 
is  characterized  by  a  vector  (p,  0,y),  where  (p,e)  are  the  polar 
coordinates  of  its  center  and  y  is  its  diameter.  The 

specification  of  disks  location  by  polar  coordinates  is  a  matter 
of  convenience  and  could  be  replaced  in  seme  application  by 
Cartesian  coordinates.  Let  5  *  {(p,  0,y) ;  (KpC®,  — iKeKrr#  0<a<y<b} 
denote  the  sample  space  of  all  possible  shadowing  disks.  Let  8 
denote  the  Borel  a  field  on  S.  a  set  B  in  5  represents 
geometrical  conditions  on  the  location  and  the  size  of  disks.  We 
consider  a  countable  collection  of  disks.  Let  N{B} ,  for  any  B 
e  8,  denote  the  timber  of  disks  having  coordinates  in  B,  i.e. , 
satisfying  geometrical  conditions  specified  by  B.  N{  * }  is  a 
0 -finite  measure  on  8.  We  consider  stochastic  structures  in  which 
N{  *}  are  random  measures.  In  these  cases,  for  each  Be8,  N{B}  is 
a  random  variable.  In  the  case  that  N{B}  has,  for  each  Be8,  a 
Poisson  distribution,  the  spatial  stochastic  process  is  called  a 
Poisson  random  field  of  shadowing  disks.  In  a  Poission  random 
field,  the  conditional  distribution  of  (N{B,}f  ...,  N{B  }),  given 
N{S}*n,  is  multinomial,  for  every  finite  partition  {Bj,  ...,  B^} 
of  S.  Accordingly,  if  the  nrnber  N{S }  in  a  Poisson  random  field 
of  disks  is  known,  we  speak  of  multinomial  fields.  In  the 
present  study  we  discuss  Poisson  random  fields. 


Hie  Poisson  random  field  model  is  obviously  an  idealization 
of  natural  phenomena.  For  example,  on  a  planar  cut  of  a  forest, 
trunks  of  trees  are  represented  by  random  disks  rather  than  by 
their  actual  random  shape.  Additional  parameters  could  represent 
other  shapes  of  shadowing  objects,  like  ellipses,  triangles,  etc. 
Furthermore,  the  Poisson  random  field  models  assumes  that,  given 
N{  S }  =  n,  the  centers  of  random  disks  are  independently 
distributed  in  the  plane.  Accordingly,  with  certain  (usually 
negligible)  probabilities  disks  may  overlap.  Clearly,  in 
practical  applications  the  fit  of  the  Poisson  model  to  the 
empirical  data  should  be  verified.  In  Poisson  random  fields  N{B} 
has  a  Poisson  distribution  with  mean 


(1.1)  u{B}  = 


■*/// 


dG(y|  p,0)H(dp,d6) , 


where  XH(dp,d6  )  is  the  expected  number  of  centers  located  in  the 
rectangle  [p,  p  +dp)  x  le,  0+dS»  G(y|  p,e)  is  the  conditional 
c.d.f .  of  the  diameter  Y,  given  the  location  of  the  center  (  p,0) . 
If  Y  is  independent  of  (p  ,0 ) ,  and  if  the  centers  are  uniformly 
distributed,  i.e.,  H(dp,de)  *  pdpde,  we  say  that  the  random  field 
is  standard.  In  this  case 


U 

(1.2)  v{B}  *  xj  H{B(y)  }dG(y) , 


where  H{B(y) }  is  the  area  of  the  subset  of  B  containing  disks 
with  a  given  diameter  y.  Standard  Poisson  and  multinomial 
fields  are  not  the  only  interesting  ones.  In  many  applications. 
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scattering  of  centers  of  disks  follow  other  distributions  than 
uniform.  In  particular,  clusters  of  disks,  each  one  following  a 
different  bivariate  normal  distribution  are  of  interest  in  various 
applications  (ecological  problems,  etc.). 
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2.  Visibility  of  Individual  Targets 

A  target  is  represented  in  the  present  paper  by  a  point  P  in 
the  plane  whose  location  is  specified  by  the  polar  coordinates 
(r,s).  The  observer  is  located  at  the  origin  0.  The  point  P  is 
visible  from  0  if  the  line  segment  OP  is  not  intersected  by  any 
random  disk.  A  random  disk  interacts  the  line  segment  OP  if  its 
coordinates  belong  to  the  set 

(2.1)  B(r,s)  «  (p,0,y);  (p/0)  e  B(r,sfy) ,  a  <  y  <  b}, 

where  B (r.s,y)  is  the  set  of  points  having  distances  fran  OP  smaller 
than  y/2  (see  Figure  1) .  A  disk  which  intersects  OP  is  said  to 
cast  shadow  on  P.  A  natural  requirement  for  shadowing  processes 
is  that  the  source  of  light  0  is  uncovered.  We  therefore 

a* 

introduce  the  structural  condition 

(2.2)  ,  C  ■  {(p,0,y);  (P/0)  e  C  (y),  a<y<b}, 

w  O 

where 

(2.3)  CQ(y)  «  {(p/0) ;  ^<p<«,  -  ir  <  0  <  tr  }. 

Figure  1  provides  a  graphical  illustration  of  the  set 
B(rrs,y)neo(y) .  A  point  P  »  (rfs)  is  visible  if  N{B(rfs)ncj  =0  . 
The  probability  that  P  is  visible  isr  therefore/ 

(2.4)  ^(r/S)  -  exp  {-v{B(r,s)ri  CQ}}. 
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In  the  standard  case,  the  area  of  B(r,s,y)Oc  (y)  is  yr.  Hence, 

b  ° 

v{B(r,s)^"|C:  }  =  Ar£,  where  £  ■  /  ydG(y)  is  the  expected  diameter 
°  a 

of  a  shadowing  disk.  In  this  case,  the  probability  that  p  is 

visible  is  (r,s)  *  exp  (-Ar  £)  irrespectively  of  the 

orientation  parameters.  A  line  of  sight,  L  ,  is  the  set  of  all 

points,  in  direction  s,  which  are  visible  from  0.  The  length  of 

L  ,  | |L  1 1 ,  has  a  c.d.f. 
s  s 


(2.5)  Pr{  | \La\ \<x]  =  l-ipjlCrS) 


Ls  is  negative  exponential 


in  the  standard  case,  the  distribution  of 
with  parameter  p=A£  (see  also  Feller  [4,p.l0] ) .  One  can  generalize  the 
above  results  to  cases  where  the  target  P  is  randomly  located.  For 
Far  example,  in  the  standard  case,  if  FR  (r)  is  the  c.d.f.  of  the  distance 
of  P  from  0,  the  probability  that  P  is  visible  is 


(2.6) 


i  =  f  e-rpdF  (r) . 


3.  Simultaneous  Visibility  of  Several  Targets 

Let  Py  ...,  Pn,  n>2,  be  arbitrary  points  in  the  plane  and 
let  0  designate  the  origin.  The  polar  coordinates  of  P^are 

(t-j,  Sj)#  j  =  If  ...,  n.  These  n  points  are  simultaneously 
visible  if  and  only  if  no  one  of  the  line  segments  (5P_.  ,  j=l, 
...»  n  ,  is  intersected  by  a  random  disk.  Let  r  ■  (r^r  *..r  t  ), 
s  *  (s  1#  ...r  sn)  and  consider  the  set  B(r,s,y)  of  points  having 
distances  not  exceeding  y/2  from  any  one  of  the  line  segments  OP_. , 
j  *  1,  n  (see  Figure  2).  Let 

(3.1)  B_ (r f s)  ■  { ( p# e,y) t  (pfflleB  (r *s,y) ,  a<y<b} . 

*1  ^  n  ~  ~  — 

The  expected  number  of  disks  which  cast  shadows  on  one  or  more  of 
the  specified  points  is  v{Bn  (r,s)nCQ  }.  Accordingly,  the 
probability  that  P  ,  ...,  Pn  are  simultaneously  visible  is 

(3.2)  >^(£,3)  s  atp{-v{Bn(r,s)n  Cq}} 

We  provide  general  formulae  for  the  determination  of  v{Bn(r,s)D  CQ} 
for  the  case  where  the  n  points  Pj  ,  j  ■  1,  n,  satisfy  the 
condition 

(3.3)  <  sx  <  s2  <  ...  <  sn  <  tt 

Condition  (3.3)  means  that  no  two  points  are  on  the  same  ray.  in 
addition  we  require  that  all  disks  satisfy  seme  condition  C, 
where  Cq  Cq,  and  v{C}  <».  Disks  which  belong  to  the  set  C  may 


or  my  not  cast  shadows  on  the  specified  points.  In  order  to 
compute  the  expected  number  of  disks  which  cast  shadows  on  at 
least  one  point  we  subtract  from  the  total  expected  number  of 
disks ,  v{C} ,  the  expected  number  of  disks  which  do  not  cast 
shadows  on  any  one  of  the  n  specified  points.  For  this  purpose 
we  define  2n  disjoint  sets  B+(ri>si,t)flc  and  B_(r.fs.,t._1)nc; 
i-1,  ...,  n,  where 


s  -s, 

"  fcn  “  * - 2  ' 


(3.4) 


t  _  3i+l~Si  *ml  _  , 

■"  2  f  X  “*  X  f  •  •  •  f  11™  X  • 


B+(r^fs^rti)D  C  is  the  set  of  disks  having  orientation 

coordinate  0e  [si  f  si+ti)  which  do  not  intersect  the  line  segment 

OP. .  The  t.  values  were  chosen  as  bisectors  of  two  adjacent 
1 

directions  in  order  to  avoid  the  possibility  that  a  disk 
belonging  to  B  +(^,  s^ft^  )H  C  will  intersect  ary  other  line 
segment  OP- ,  j^L.  The  sets  B_(r. ,  s.  ,t.  JO  c  have  similar 

1  1  1*“X 

interpretation ,  for  disks  in  C  with  orientation  coordinates 
0e[s^-t^_^#  s^) .  The  sets  B±(r,s,t)  of  disks  parameters  (p,  8, y)  are 
described  in  terms  of  sets  of  center  ponits  ( p,  0)  of  disks  having 
a  given  diameter  y.  These  sets  are  denoted  by  B±  (r,s,t,y) .  In 
Figure  2  we  illustrate  a  partition  of  the  ( p,  6)  plane  in  terms  of 
such  sets  for  a  given  y.  Accordingly 


(3.5)  B±  (r,s,t)-{(p,0,y);  (p»0)  eB±  (rfs,t,y) ,  a<y<b}. 


Figure  3  provides  geometrical  details  for  the  definition  of 
B±(r,s,t,y)  .For  a  formal  definition,  let  t(r,y)  *  min(t,tan’"1  (y/2r) ) 

and  let 


r  +  y/2 


,  if  4>  =  0 


(3.6)  P(r,<|>,y)  *  <  r[cos  <t>  +  ( (y/2r) 2-sin2  )1/2],  0«Kt(r,y) 


y/  (2sin<t>) 


,  t(r,y)«jKt. 


Then, 

(3.7)  B+(r,s,t,y)  -  {(P,9);  P(r,0-s,y)<p«»  and  0  <9-s<t}. 
Similarly, 

(3.8)  B_(r,s,t,y)  ■  {(p,0);  p(r,s-6,y)  <  p<~,  and  O<s-0<t>. 
The  expected  ranter  of  disks  in  B±(r,s,t)n  CQ  is  denoted  by 

A  K±  (r,s,t).  Me  formally  define  K±  (r,s,t)  ■  0  for  all  UO. 

Generally,  for  t>0, 

(3.9)  K±(r,a,t)  -  JJj  dG(y|p, 6>H(dp,d9) . 

B±(r,s,t)DC 


The  above  K-f  unctions  can  be  used  to  express  the  probabilities 
<>n(r,s)  more  explicitly.  The  expected  matter  of  disks  in  C,  which 
oast  ^radons  an  one  or  sore  of  the  n  specified  points,  is 


(3.10)  v{Bn(r,s)ncMC}-xij^[KJri,si,ti_1)+  VWV1' 
Substitution  of  (3.10)  in  (3.2)  yields  more  explicit  formula  for 


SBSgm 


Consider  a  target  moving  at  a  constant  speed  along  a  curve  C 
in  the  plane.  We  are  interested  in  determining  the  distribution 
function  of  the  proportion  of  time  in  which  the  target  is  not 
shadowed  by  random  obs  cur  ring  elements  while  moving  between 
points  PL  and  Pu  on  C.  This  proportion  of  time  is  equivalent  to 
the  proportion  of  the  length  of  the  visible  segments  of  C  between 
Pr  and  P  .  We  assume  that  C  is  a  star-shaped  curve.  Such  a  curve 
intersects  any  ray  from  the  origin  at  most  once.  Hie  curve  C  is 
specified  by  a  function  r(6)/  s.j<9_<su,  which  yields  the  distance 

of  C  from  0  in  direction  6.  We  further  assume  that  r(s)  is 
piecewise  continuously  differentiable .  Let  sL  and  su  be  the 
orientation  coordinates  of  PT  and  P  .  respectively.  The  length 

Ml 

of  C ,  between  these  two  points'  is 


l(s)ds  / 


(4.1) 


where 


(4.2)  £<s)  =  [r2(s)  +  (gf  r(s))2]1/2 


The  proportion  of  the  visible  segments  of  C,  between  and  P^, 


(4.3) 


i  A 


I(s)£(s)ds  , 


where  I (s)  =  1  if  the  point  (r(s),s)  on  C  is  visible,  and  I(s)  =0 
otherwise.  V  is  a  random  variable  assuming  values  in  [01]. 
Hie  distribution  of  V  is  a  mixture  of  a  two-point  distribution 


concentrated  on  {0,1}  and  a  distribution  F^*(v)  on  (0,1).  Under 
some  general  conditions  the  distribution  concentrated  on  (0,1)  is 
absolutely  continuous.  It  sufficies,  for  example,  that  G(y|p,6) 
is  absolutely  continuous,  or  that  H(dpvd6)  =  h(p ,  6)dpde, 
h(P,  9)_>0.  Let  pQ  *  Pr{V*0}  and  Pl  *  Pr{V*l}.  The  c.d.f  V  of 
can  be  generally  represented  by  the  formula. 

0  ,  v  <  0 

(4.4)  Fy(v)  =  PQ  +  (l-p^-pj)  F£  (v)  ,  0  £  v  <  1 

1  ,  1  £  V  . 

The  probability  that  the  section  of  C,  between  PT  and  P„  ,  is 
completely  visible  is  given  by 

(4.5)  Pl=exp{-[v{C}-A[K+(r(su),  su,t)+K_(r  (sL)  ^t)  ]  ]  } 


where  t  =  tt -(su~sl)/2  and  C  is  the  set  of  shadowing  disks. 

AK+(r(sy  ),  su,t)  is  the  expected  number  of  random  disks  having 
orientation  coordinate  0e(s  t  ,s +t) ,  which  do  not  cast  shadows  on 
the  segment  of  C  under  consideration.  Similarly,  AK_(r(sL)  rfyt) 
is  the  expected  number  of  random  disks  with  orientation 
coordinate  9e(s  j^-t^  ),  which  do  not  cast  shadows  on  the 

intersecting  segment  of  C .  Explicit  expressions  for  the 

probability  of  complete  invisibility  pQ  and  of  the  c.d.f.  Fv(v) 
are  not  available.  An  approximation  to  F  (v)  by  a  beta 


distribution  will  be  discussed  in  Section  4.2.  A  discrete 
approximation  will  be  studied  elsewhere. 


4.1  Moments  of  V 

The  n-th  moment  of  V  is,  according  to  (4.3) , 


(4.6)  u  =  —  E  {  (  /  6(s)I(s)ds)n) 

L  *L 

SU  ,su  n  n 

...  f  E{tt  I(s. )}  tt  £(s. )ds. 

J  i=l  1  i-1  11 


= 

Ln  J 

-/ 
A  J 

n 

E{tt  I(s.  ) 
i=l  1 

n 

}  TT  £(s.  )ds.  . 
i=l  1  1 

= 

SL  f 

Ln  ' 

-/ 
A  J 

Pn(31,  •• 

n 

.,  s  )  tt  l{s.)ds.  , 
n  i=l  1  1 

where  A  = 

■ 

(0 

M 

A 

• 

• 

£  s^}  and 

• 

n 

(4.7) 

pn(s 

^  •••» 

»_>  - 

n 

*n(5'5}  * 

E{  tt  I(s^)}. 

Applying  the  K-f unctions,  which  are  defined  in  (3.7) ,  one  can 
express  the  probability  of  simultaneous  visibility  E^(s  lf  ...,  sn) 
in  the  form 

(4.8)  pn(slf  aj  -  exp{^{C}}. 
n 

.expU  E  (K  (r(s.),s. ,t.  .)  +  K  (r. (s. ) ,s. ,t. ) ] } 

A  AAA  T  1  A  A 

n-1 

"  exp{-v{C>}A0(s.,sn)  JT  A(s.,si+1), 


(4.9)  Ao(x,y)=exp{A[K_(r(x),x,Tr-  J^-)+K+(r  (y)  ,y,n  -  *^)]} 


-7T  £  x  f_  y  _<  tt,  and 


(4.10)  A(x,y)*exp{A[K+(r(x),x,^)  +  K_(r(y)  ,y,^)]  },  -it  £  x  £  y  <tt. 


Substitution  o£  (4.8)  and  (4.6)  yields  a  recursive  formula  for 
the  moments  of  V.  Let  £*(s)  a  £  (s)/L  and  define  recursively. 


(4.11) 


Qjbc^y)  =  A(x,y) , 


Qi+l(3C,y)  =  (i+l)^f  £*(z)A(z,y)Qi(x,z)dz  ,  i*l,2,, 


Finally,  the  moments  of  V  are 

,su 

Uj^expt-viC}}  J  £*(x)AQ(x,x)dx  , 

(4.12)  and  for  n>2 

s  s 

ru  ru 

Vn  *  n  exp{-v{C>>  J  £*(x)  J  £*(y)AQ(x/y)(^_1(x,y)dy  dx. 
®L 

Notice  that  if  su  -sL<  tt  then  AQ(x,y)  ■  A^  (x)A2(y)  and  expression 
(4.12)  assumes  a  simpler  form.  Formulae  (4.11)  and  (4.12)  are 
convenient  for  the  recursive  determination  of  the  moments. 
Finally  we  remark  that  the  sequence  {  un  ;n>l}  is  monotonically 
decreasing  and  lim  un  »  p ,  .  Indeed,  according  to  the  dominated 

Tt*» 

convergence  theorem. 


(4.13) 


Un  ^  -  Unfll^-pj)  /  und^(u)4p  ] 

rr*»  n-*»  0 

1  ° 

“  (l-P0“Pi)  f  (lim  un)dF*  (u)+p. 

J  rv4CO  V  A 
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4.2  A  Beta-Mixture  Approximation  to  the  CDF  of  N. 

In  order  to  approximate  the  c.d.f.  Fv(v)  we  consider  a 
mixture  of  a  two-point  distribution  and  a  beta-distribution. 
This  approximation  is  given  by 

t 

0  ,  v  <  0 

*0 

(4.14)  Fv(v)W  pQ  +  (l-po-p1)  Iv(a,B)  ,  0  <  v<l 

„ 1  >  1  <  v  , 

X 

where  Ix(a,b)  =  ^  J  ua  ^(l-uj^du  is  the  incomplete  beta 

o 

function  ratio,  0<a,  b<  00 .  The  parameters  pQ  ,  a  and  B  are. 

determined  by  equating  the  first  three  moments  of  V  to  those  of 

*** 

Fv(V).  The  moment  equations  yield  the  formulae 

PQ  =  1-Pjl  -  (Uj^-Pi)  Cot+B)  /a  , 

(4.15)  a  =  (2(y*)2-y*(y*+y*))/D 
B  *  (yj-y*) (y£-y£)/D, 

vtere  y*  -  and  D  -  y*  y*  -  (y*)2. 

This  approximation  was  previously  applied  by  Yadin  and  Zacks  [12] 
in  a  related  shadowing  problem.  In  Section  6  we  provide  a 
numerical  example  which  illustrates  this  approximation  in  a 
special  case. 


The  Evaluation  of  the  K-Functions  Under  a 


ial  Field 


Spec 

Structure 

The  evaluation  of  the  R-f unctions  depends  on  the  particular 

structure  of  the  field  of  random  disks.  In  the  present  section  we 

consider  a  class  of  field  structures,  which  is  prevalent  in  many 

applications,  and  Which  yields  relative  simple  formulae.  This 

class  is  characterized  by  the  condition  that  centers  of  disks  are 

scattered  within  a  regions  which  is  situated  between  0  and  the 

observed  points  P,  . . . ,  Pn,  or  the  curve  C . 

The  regions  in  which  centers  of  disks  are  scattered  are 

specified  by  two  star-shaped  curves  U  and  W  between  the  origin  and 

the  target  curve  (or  points) .  These  curves  are  specified  by  two 

star-shaped  functions  u(9)  ,w(0)  ,0e[0_  #0  ],  respectively,  and  the 

L  U 

regions  is  given  by 

(5.1)  S  =  {(p,0);  u(0) <p<w(0) ,  9_  <  0  <  0  }  . 

—  i-i  ““  U 

The  requirement  that  disks  having  centers  in  the  region  S  and 
diameter  in  [a,b] ,  would  not  cover  the  origin  or  intersect  target 
points  implies  the  following  conditions  on  u(6)  and  w(e ) : 

(5.2)  |<u(0)  ,  0e[0L,0u]'  ; 

and 
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Fig.  4.  The  Geometry  Of  Disks  Under  A  Special  Field  Structure 


(5.3) 


w(9)  oos(9-s)  +  ((^)2  -  w2(9)sin2(9-s)]Js<r(s) , 

for  all  9e[9L/9ul , 

where  (r(s),s)  is  any  target  point.  Condition  (5.3)  insures  that 

the  distance  r(s)  between  0  and  (r (s) ,s)  is  greater  than  the 

distance  £ (s)  between  0  and  P  =  (S (s) ,s) ,  where  P  is  a  point  on  a 
\  ^ 
circumference  of  a  disk  of  the  maximal  diameter,  b,  whose  center 

(w (9), 9)  is  on  (I)  (see  Figure  4).  The  structural  condition  imposed 

on  the  field  of  random  disk  is  C  *  S  x  [a,b] .  In  the  following 

development  we  assure,  for  the  sake  of  simplifying  notation  that 

the  centers  of  disks  are  distributed  within  S  according  to  an 

absolutely  continuous  measure,  i.e. , 


(5.4) 


H(dp,d9)  = 


h(p,9)dpd9 

0 


,  if  (p,9)eS 


,  otherwise. 


A  disk  centered  at  (P  ,9)  does  not  intersect  a  ray  from  0  having 
orientation  s,  if  its  diameter  y  is  smaller  than  2psin| 9-s| ,  or  if 
1 9  -s|  >  tt/2.  Let  G*  (p  ,9  ,s)  denote  the  probability  of  this  event, 
i.e. 


G(2psin|9-s||p,9)  ,  if  |9-s|<sin  H 


(5.5)  G*(p,9,s) 


,  otherwise 


Let 


w(9) 


(5.6)  k(s,9)  =  J  G* (p,9,s)  h(p,9)dp  . 


■  ^ 


6.  An  Example 


In  the  present  example  we  develop  all  the  previous  formulae 
explicitly  and  provide  numerical  illustrations  for  the  following 
simple  case. 

(1)  Hie  Poisson  field  is  standard,  with  intensity  x; 

(ii)  Hie  diameters  of  disks  are  uniformly  distributed  in 
[a,b],  0<a<b; 

(iii)  Hie  target  curve  C  is  an  arc  on  a  circle  of  radius 
r,  specified  by 

C  =  { (p,0) ;  p=r  ,  <  0  <  suJ  . 

(iv)  Hie  region  S  is  an  annular  region  bounded 

between  two  concentric  arcs,  U  and  to  of  radii  u 

and  w,  respectively,  which  satisfy  the  condition 

b  <  u  <  w  <  r-b/2. 

2 

That  is 

S  »  { (p, 0)  ?  u  <  p  <  w;  0_  <  0  <  0,,} 
and  t-I-fl- 

Hie  ejected  number  of  disks  in  S  obtains  the  siitfpie  expression 


0  w 

U 

0LU 


pdpd0  =  -|(w2-u2)  (0u-0L) 


(6.1) 


v{C}  *  X 


To  develop  the  K-functions  we  substitute  in  (5.6)  the  function 


(6.2) 


,  if  1 6-s  |  <sin-1  (-^) 


G*(p,s,0)  =  < 


2psin  0-s  -a  a*  -1,  b» 

- bia"1  ,  if  sin  (^)<|0-si<sin  (^) 


2p 

,  if  sin-1  fe^)  <_  1 0-s  | . 


2p 


Hius  one  obtains 


(6.3)  k(s,0)  =  kw(ls-0|)-lcli(|s-0|)  , 

where,  for  t>0 


(6.4)  kv(t) 


[4  v3sxnt-  |  v2+  — a  -  ]  /  (b-a) 

J  *  24sint 


I.  2  a2-fab+b2  , 

2 lv  “  2  1 
*  12  sin  t 

■j[v2-  (a2+ab+b2)  /12] 


,  if  t  sin-1  (~) 

,  if  sin-1  (^)£t£sin_1  (-^] 


,  if  sin-1  (-^)  <  t  < 


,  if  tt/2  <  t. 

For  positive  t  values  not  exceeding  minte-s,  s-e  ),  the 

u  L 

K-functions  satisfy 


ro|=t 


(6.5)  K  (r,s,t)  *K  (r,s,t)  -  K*(t)  ■  K*(t)  -K*(t), 
T  —  w  u 


where 


(6.6) 


<(t)  =  I  k 
v  j  v 


(T)dT 


0 

,  if  t  <  sin 

K(1)  (t) 

,  sin”1  (J)  £  1 

\(2)  (t) 

'  sln~1 (2v}  i 

kv<3)  (t) 

,  tt/2  £  t 

-I,  b. 
1  2v 


in  which 


(6.7) 


*v]  <«  *  (x  f*v2-a2)1/2  -  2v  cos(t)]  - 


2 

f  V2  [t-sin-1  (^) ]  +  Jj.  [4v2-a2)1/2  -  a  cotan(t)]}. 


(6.8) 


-K^W1^))  +f  (t-sirf1^)) 


2v"  T  2 
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(6.9) 


(3)  ,  v  ‘21  (I,  +  (t-  §>  (£  -  2^*!  ,. 


V"  <«  -  \ 
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If  ...,  Pn  are  points  on  C  such  that  sL<s1<s2<  ...  <  sn<su 
then,  the  probability  of  simultaneous  visibility  is  given  by  the 
formula 


(6.10)  Pn(sl'  ***/  sn^  =  (-  |-(w2-u2)  (eu-eL) ) 

n-1  s-.-s. 

.expU[K*(s1-0L)  +  K*(8u-sn)  +  2  Z  K*(-|-— ■ ■■)]}. 
Furthermore,  the  probability  that  C  is  completely  visible  is 

(6.11)  p^  *  exp{-  ^(w2-u2)  (0U-0L) }. 
exp(X[K*(sL-0L)  +K*(0u-su)] 


In  liable  6.1  we  provide  numerical  values  of  the  probabilities 

p  (s, ,...,  a  )  and  p.  for  the  following  case.  The  region  s  is 
n  1  n  1 

located  in  an  annular  region,  with  e_  ■  -  it/2,  0  =  ir/2,  u  -  .5, 

Xj  u 

w  *  .75.  The  distributions  of  Y  is  uniform  on  the  interval  (.1, 
.3).  The  intensity  of  the  Poisson  field  is  X  -  1(2)9.  n  target 
points  are  placed  on  a  half-circle  with  radius  r  ■  1,  and 
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U’a  ‘^  '4 


.34786 


.25726 


Table  6.1.  Probabilities  of  Simultaneous  Visibility  of  n  Points 


Uniformly  Placed  on  a  Half  Circle;  u  *  .5,  w  *  .75,  r  *  .1,  b  * 


For  tiie  computation  of  moments  of  V,  we  consider  a  target  curve  C , 
on  the  half-circle  of  radius  r  *  1  with  s^5*  -tt/3  and  su  =  i/3.  In 
this  case,  l*(z)  for  -^<_z  <_v/3,  AQ(x,y)  *  A1(x)A2(y), 

where  A^x)  ■  exp  {XK*(x  -tj)},  and  ^  (y)  -  exp  (XK*f£-y)}, 

- 5- *  <  y  - tt/3  •  sijailarly'  A(x,y)  -  exp(2XK* £=£)}. 

- 1  <  x  <  y  <  ir/3.  Since  A  (x,y)  is  factored  to  the  product  of 


(x)A2(y) ,  the  recursive  integration  needed  for  the  moments 
obtains  the  following  form. 


Let  Qq(x)  =  A1(x) 

and  for  j  >  1,  let 

(6.13)  Q.(x) 

The  moments  of  V  are  determined  by  the  formula 


=  -§j  J  Qj_i  (y)A(y,x)dy  . 


(6.14)  yn=exp{-.15625X7r}  J  (^(yjAty^cty. 

*L 

The  probability  that  C  is  completely  visible  is  given  by 

(6.15)  P1=«q?{- .  15625Xit  }exp{  2  AK*  (ir/6 ) ) . 


The  integrations  in  (6.13)  and  (6.14)  were  performed  numerically 
In  Table  6.2  we  present  the  values  of  un  for  n  ■  1(1)10  and  the 
value  of  p1  for  the  same  parameters  as  in  Table  6.1. 


1 

.951 

.912 

.880 

.854 

.832 

.814 

.799 

.787 

.776 

.767 

.686 

3 

.860 

.758 

.682 

.623 

.577 

.541 

.511 

.487 

.467 

.451 

.322 

5 

.778 

.631 

.529 

.456 

.402 

.360 

.328 

.303 

.282 

.266 

.152 

7 

.703 

.526 

.412 

.335 

.280 

.241 

.212 

.189 

.171 

.157 

.071 

9 

.636 

.438 

.321 

.246 

.196 

.162 

.137 

.118 

.104 

.094 

.033 

Table  6.2  Moments  of  V  for  X=l(2)9  and  ^  =  -  tt/3,  =  ir/3,  0L  -  it, 

9u  =  Tr/2,  u  =  *5/  w»  .75,  r  =  1.0,  a  =  .1,  b  =  .3 


In  Table  6.3  we  present  the  moments  of  V  and  the  corresponding  manents 
of  the  mixed-beta  approodmation,  for  a  slightly  different  case,  in 
which  C  is  an  arc  an  a  circle  of  radius  r  =  1,  between 
a,.  =  -tt/18  to  s  =  -s, . 


Intensity 


raxiroation  (lower)  for  the  case  of  sL 
.5,  w  =  75,  r  *  1.0,  a  =  .1,  b  =  .3 


Notice  that  is  the  same  in  tables  6.2  and  6.3,  for  all  A 
values.  The  moments  of  the  mixed-beta  approximation  are  generally 
close  to  those  of  V.  this  shews  an  adequate  approximation.  Hie 
values  of  vn,  n  >  1  should  be  monotonically  decreasing.  Hie 
slight  deviation  for  high  order  moments  and  small  A  is  due  to 
numerical  errors  in  the  integration.  In  Table  6.4  we  present  the 
parameters  of  the  mixed-beta  approximating  distribution 
corresponding  to  Table  6.3. 


.1733 

.0020 

.9005 

1.1161 

1.0398 

.2765 

.0152 

.7302 

1.1405 

1.0967 

.3290 

.0326 

.5921 

1.0956 

1.1194 

.3590 

.0536 

.4801 

1.0450 

1.1386 

.3754 

.0768 

.3893 

0.9937 

1.0967 
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